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Abstract 

Wc introduce a discrete-time fractional calculus of variations on the time scale 
/iZ, /i > 0. First and second order necessary optimality conditions are estab- 
lished. Examples illustrating the use of the new Euler-Lagrange and Legendre 
type conditions are given. They show that solutions to the considered frac- 
(~| , tional problems become the classical discrete-time solutions when the fractional 

"*li ' order of the discrete-derivatives are integer values, and that they converge to the 

fractional continuous-time solutions when h tends to zero. Our Legendre type 
condition is useful to eliminate false candidates identified via the Euler-Lagrange 
fractional equation. 
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1. Introduction 



The Fractional Calculus (calculus with derivatives of arbitrary order) is an 

^^ , important research field in several different areas such as physics (including 

\^ ' classical and quantum mechanics as well as thermodynamics), chemistry, biol- 

. 5t 1 ogy, economics, and control theory [3|, |lO|, |40|, |42, |48[. It has its origin more 



than 300 years ago when L'Hopital asked Leibniz what should be the meaning 
of a derivative of non-integer order. After that episode several more famous 
mathematicians contributed to the development of Fractional Calculus: Abel, 
Fourier, Liouville, Riemann, Riesz, just to mention a few names [30|,|47[. In the 
last decades, considerable research has been done in fractional calculus. This is 
particularly true in the area of the calculus of variations, which is being subject 
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to intense investigations during the last few years [lOJ, |ll|, |4J, |45|. Applica- 
tions include fractional variational principles in mechanics and physics, quan- 
tization, control theory, and description of conservative, nonconservative, and 
constrained systems 10|, ll5|, |l6|, |45| . Roughly speaking, the classical calculus of 



variations and optimal control is extended by substituting the usual derivatives 
of integer order by different kinds of fractional (non-integer) derivatives. It is 
important to note that the passage from the integer/classical differential calcu- 
lus to the fractional one is not unique because we have at our disposal different 
notions of fractional derivatives. This is, as argued in [10|, |4J|, an interesting 
and advantage feature of the area. Most part of investigations in the fractional 
variational calculus are based on the replacement of the classical derivatives 
by fractional derivatives in the sense of Riemann-Liouville, Caputo, Riesz, and 
Jumarie [l|, y, |lO|, |27|. Independently of the chosen fractional derivatives, one 



obtains, when the fractional order of differentiation tends to an integer order, 
the usual problems and results of the calculus of variations. Although the frac- 
tional Euler-Lagrange equations are obtained in a similar manner as in the 
standard variational calculus [44|, some classical results are extremely difficult 
to be proved in a fractional context. This explains, for example, why a fractional 
Legendre type condition is absent from the literature of fractional variational 
calculus. In this work we give a first result in this direction (cf. Theorem 13. 6p . 
Despite its importance in applications, less is known for discrete-time frac- 



tional systems [4J]. In [39| Miller and Ross define a fractional sum of order 
u > via the solution of a linear difference equation. They introduce it as (see 
^for the notations used here) 

^-fit) - w^ E(^ - <'^)Y''-'^fis)- (1) 

^ ^ s—a 

Definition ([T]) is analogous to the Riemann-Liouville fractional integral 

of order v > 0, which can be obtained via the solution of a linear differential 
equation 39l|40|. Basic properties of the operator A^'^ in ([T|) were obtained 



in [39|. More recently, Atici and Eloe introduced the fractional difference of 
order a > by A"/(i) = A™(A"~™/(i)), where m is the integer part of a, 
and developed some of its properties that allow to obtain solutions of certain 
fractional difference equations [8|, |9| . 

The fractional differential calculus has been widely developed in the past few 
decades due mainly to its demonstrated applications in various fields of science 
and engineering 30|, |40|, |43[ . The study of necessary optimality conditions for 



fractional problems of the calculus of variations and optimal control is a fairly 
recent issue attracting an increasing attention - see [ll l^; LE l22|, |23|, |25|, |26|, |41 1 
and references therein - but available results address only the continuous-time 
case. It is well known that discrete analogues of differential equations can be 
very useful in applications 13, _3l|, |29[ and that fractional Euler-Lagrange dif- 
ferential equations are extremely difficult to solve, being necessary to discretize 
them [^, |ll|. Therefore, it is pertinent to develop a fractional discrete-time 
theory of the calculus of variations for the time scale {h'E)a, h > (cf. defini- 



tions in Section [2]). Computer simulations show that this time scale is particu- 
larly interesting because when h tends to zero one recovers previous fractional 
continuous-time results. 

Our objective is two-fold. On one hand we proceed to develop the theory of 
fractional difference calculus, namely, we introduce the concept of left and right 
fractional sum/difference (cf. Definition 12. 8p . On the other hand, we believe 
that the present work will potentiate research not only in the fractional calculus 
of variations but also in solving fractional difference equations, specifically, frac- 
tional equations in which left and right fractional differences appear. Because 
the theory of fractional difference calculus is still in its infancy [8|, |9|, |39| , the 
paper is self contained. In Sj2]we introduce notations, we give necessary defini- 
tions, and prove some preliminary results needed in the sequel. Main results of 
the paper appear in SJ31 we prove a fractional formula of /i-summation by parts 
(Theorem 13. 2p , and necessary optimality conditions of first and second order 
(Theorems 13.51 and 13.61 respectively) for the proposed /i-fractional problem of 
the calculus of variations (jTTl) . Section |4] gives some illustrative examples, and 
we end the paper with ^Dof conclusions and future perspectives. 

The results of the paper are formulated using standard notations of the 



theory of time scales [20|, [32, l33| . It remains an interesting open question how to 
generalize the present results to an arbitrary time scale T. This is a difficult and 
challenging problem since our proofs deeply rely on the fact that in T = (hl^a 
the graininess function is a constant. 



2. Preliminaries 

We begin by recalling the main definitions and properties of time scales 
(cf. [ISJ, 1201 and references therein) . A nonempty closed subset of R is called 
a time scale and is denoted by T. The forward jump operator cr ; T — > T is 
defined by a{t) = inf {s G T : s > i} for all t G T, while the backward jump 
operator p : T — !> T is defined by p{t) — sup {s G T : s < i} for all i G T, with 
inf = supT (i.e., cf{M) ~ M if T has a maximum M) and sup0 — inf T 
(i.e., p(rn) = m if T has a minimum m). A point t G T is called right-dense, 
right- scattered, left-dense, or left- scattered, if a{t) — t, a{t) > t, p{t) = t, or 
p{t) < t, respectively. Throughout the text we let T = [a, &] n T with a < b and 
T a time scale. We define T** ~ T\{p{b),b], T'^ — (T'^)'" and more generally 
!"*■ = IT'' J , for n G N. The following standard notation is used for a 

(and p): a^{t) ~ t, cr^^{t) — [a o cr"^^)(<), n G N. The graininess function 
M : T -> [0, oo) is defined by p{t) = cr(t) - 1 for aU t G T. 

A function / : T — > R is said to be delta differentiable at t G T** if there is 
a number f^{t) such that for all e > there exists a neighborhood U oi t (i.e., 
U ^ {t- 5,t + 5)r\T for some S > 0) such that 

\f{a{t)) - f{s) - /^(t)(a(t) - s)\ < e\<j{t) - s\, for aU s G U. 

We call f^{t) the delta derivative of / at t. The r^'^ — delta derivative (r G N) 
of / is defined to be the function /^"^ : T**' -^ R, provided /^' is delta 
differentiable on T** . For delta differentiable / and g and for an arbitrary 
time scale T the next formulas hold: /'^(t) = f{t) -\- ^{t)f^{t) and 

(/5)^(i) = f^it)9''{t) + f{t)9^{t) = f^{t)9{t) + f^t)g^{t), (2) 



where we abbreviate f o a hy f^ . A function / : T — ;■ R is called rd- continuous 
if it is continuous at right-dense points and if its left-sided limit exists at left- 
dense points. The set of all rd-continuous functions is denoted by Crd and the 
set of all delta differentiable functions with rd-continuous derivative by Cl^. 
It is known that rd-continuous functions possess an antiderivative, i.e., there 
exists a function F G Cj?j with F^ = /. The delta integral is then defined by 

JI^ fit)At = F{b) - F{a). It satisfies the equality /j"^*^ /(t)At = n{t)f{t). We 
make use of the following properties of the delta integral: 

Lemma 2.1. (cf. ]2^, Theorem 1.77]) Ifa,beT and f,g eCrd, then 

1. J^ficT{t))g^it)At= ifgmf=':~f^f^it)git)At; 

2. j;;f{t)g^{t)At= {fg){t)\lZl-f;:f^{t)g{a{t))At. 

One way to approach the Riemann-Liouville fractional calculus is through 
the theory of linear differential equations [43[. Miller and Ross [39| use an 
analogous methodology to introduce fractional discrete operators for the case 
T = Za = {a, a 4- 1, a + 2, . . .}, a e R. Here we go a step further: we use the 
theory of time scales in order to introduce fractional discrete operators to the 
more general case T = {hZ)a — {a, a + h,a + 2h, . . .}, a G M, /i > 0. 

For n G No and rd-continuous functions pi : T ^- M, 1 < i < n, let us 
consider the nth order linear dynamic equation 

n 

Ly = , where Ly = y^" + ^fty^" ' ■ (3) 

A function y : T ^- M is said to be a solution of equation ([3]) on T provided y is 
n times delta differentiable on T" and satisfies Ly{t) — for all t G T** . 

Lemma 2.2. fM, P- 239] // 2 = (zi, . . . , z„) : T -s- R" satisfies for all t G T'' 



z^ = A{t)z{t), where A 
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then y — zi is a solution of equation ([3]). Conversely, if y solves ([3]) on T, then 
z = (y,y^, . . .,y^"'') : T ^ R satisfies g]) for all t G T''" 

Definition 2.3. ]^2C, p. 239] We say that equation (0] is regressive provided 
I + i-i(t)A(t) is invertible for all t G T", where A is the matrix in ([4]). 

Definition 2.4. ^, p. 250] We define the Cauchy function y : T x T''" -^ R 
for the linear dynamic equation Q) to be, for each fixed s € T'^ , the solution 
of the initial value problem 

Ly^O, y^*((a(s),s)=0, < i < n - 2, y^"~\ia{s),s) = 1 . (5) 



Theorem 2.5. J^d . p. 251] Suppose {2/i,...,y„} is a fundamental system of 
the regressive equation Q). Let f 6 Crd- Then the solution of the initial value 
problem 

Ly = f{t), y^\to) = 0, 0<i<n-l, 

is given by y{t) = J y{t, s)f{s)As, where y(t, s) is the Cauchy function for ^. 

It is known that y{t, s) := i/„_i(i, a{s)) is the Cauchy function for y — 0, 
where Hn-i is a time scale generahzed polynomial [20, Example 5.115]. The 
generalized polynomials Hk are the functions H^ : T^ ^ R, fc 6 Nq, defined 
recursively as follows: 

i/o(i,s) = l, Hk+i{t,s)^ f Hk{T,s)AT, fc = l,2,... 

for all s,t E T. If we let H^{t,s) denote, for each fixed s, the derivative of 
Hk{t,s) with respect to i, then (cf. [20|, p. 38]) 

H^it, s) = Hk-i{t, s) for fc e N, teV". 

From now on we restrict ourselves to the time scale T = {hZ)a, h > 0, for 
which the graininess function is the constant h. Our main goal is to propose and 
develop a discrete-time fractional variational theory in T — {hZ)a- We borrow 
the notations from the recent calculus of variations on time scales 18|, |2J, |33| • 



How to generalize our results to an arbitrary time scale T, with the graininess 
function /i depending on time, is not clear and remains a challenging question. 
Let a £ M and /i > 0, {h'Z)a — {a, a + h,a + 2h, . . .}, and b = a + kh for some 
k E N. We have a{t) = t + h, p{t) = t — h, fi{t) = h, and we will frequently 
write fit) = f{a{t)). We put T = [a, b] n {hZ)a, so that T*" = [a, p(b)] n {hZ)a 

and T" = [a, p^(6)] n {hZ)a. The delta derivative coincides in this case with 

f^U) — fit) 

the forward /i-difference: f (t) = --r . li h = 1, then we have the 

fi{t) 

usual discrete forward difference A/(f). The delta integral gives the ft,-sum 
(or /i-integral) of /: / f{t)At = \J f{kh)h. If we have a function / of two 

J a I a 

'—h 
variables, f(t,s), its partial forward /i-differences will be denoted by At^h and 

Asji, respectively. We will make use of the standard conventions X]t=c /(*) = 



r-l 



0, c G Z, and Ili^o /(*) ~ 1- Often, left fractional delta integration (resp., 
right fractional delta integration) of order j/ > is denoted by aA^'^ f{t) (resp. 
tA^^ f{t)). Here, similarly as in Ross et. al. [46], where the authors omit the 
subscript t on the operator (the operator itself cannot depend on t), we write 
aA-V(t) (resp. hA-'^fit)). 

Before giving an explicit formula for the generalized polynomials Hk on hZ 
we introduce the following definition: 

Definition 2.6. For arbitrary x,y eM the h-factorial function is defined by 

,iy) ._ J,, ni + i) 

^ ■" ni + i-y)' 

where T is the well-known Euler gamma function, and we use the convention 
that division at a pole yields zero. 



Remark 2.1. For h = 1, and in accordance with the previous literature ^, we 
write x^y' to denote Xy^ . 

Proposition 2.1. For the time-scale T = {h1)a one has 

(t-s)^''^ 

HkU, s) -.^ -^ ^^^^ for all s,t eT and k e No ■ (6) 

fc! 

To prove (O we use the following technical lemma. Throughout the text the 
basic property r(a:: + 1) = xT{x) of the gamma function will be frequently used. 

Lemma 2.7. Let s E T. Then, for all t eT'^ one has 

i it-^tin _ it-^ 

^'^'\ (fc + l)! /-^fc!~- 
Proof. The equality follows by direct computations: 

r (t-^£^i 1 r (ow-^r^ {t-s)r'' \ 

*'''! (fc + l)! J h\ (fc + l)! (fc + l)! j 

h''+^ r T{{t + h-s)/h + l) T{{t-s)/h + l) 



h{k + 1)! i r((t + h- s)/h + 1 - (fc + 1)) r((i - s)/h + 1 - (fc + 1)) 
/i'= ({{t-s)/h + l)Tiit~s)/h+l) T{{t-s)/h+l) 



(fc + 1)! t ((* - s)/h - k)T{{t - s)/h - fc) T{{t - s)/h - fc) 
h'^ ( ri{t-s)/h+l)\ {t-s)i'^ 



fc! {T{{t-s)/h+l-k)j fc! 

D 
Proof, (of Proposition 12. ip We proceed by mathematical induction. For fc = 

^o(*' ^) - ^'^ r(^ + i-o) - r(^ + i) - ' ■ 

Assume that (O holds for fc replaced by in. Then by Lemma [^771 

i/™+i(t,s) = / iJ„,(r,s)Ar= / ^ " Ar=^^^^ 



i! (m + 1)! ' 

which is ([6]) with fc replaced by ttt, + 1. D 

Let yi{t), . . . , y„(i) be n linearly independent solutions of the linear homoge- 
neous dynamic equation y — 0. From Theorem l2.5l we know that the solution 
of © (with L = A" and io = a) is 

yit) = A-fit)= ^-^^ f{s)As^—- Y.{t-aikh)t-'^fikh)h. 

•^^ ^ ' ^ ^ k=a/h 



Since y^{a) — 0, z — 0, . . . , n — 1, then we can write that 

t/h — n 
^ ^ k^a/h /'J^ 



Note that function t -^ (A^"/)(f) is defined ioi t — a + nh mod(/i) while 
function t — )■ f{t) is defined for i = a niod(ft.). Extending ([7]) to any positive 
real value v, and having as an analogy the continuous left and right fractional 
derivatives [40[, we define the left fractional h-suui and the right fractional h- 
sum as follows. We denote by J-r the set of all real valued functions defined on 
a given time scale T. 

Definition 2.8. Let a e R, h > 0, b = a + kh with k e N, and put T = 
[a,b] n {h'Z)a. Consider f G J-j. The left and right fractional h-sum of order 
ly > are, respectively, the operators a^h" ■ -^t -^ J^f+ and ^A^" : JFj -^ J^f-, 
T^ — {t zL i/h : t € T} , defined by 

t 

,A,r/(t) ^ -i- r * '^'^\t^a{s)ir'^f{s)As^-^Y,{t^a{kh)t-'^f{kh)h 






'^^^^ = T{y)Jt+J'"'^^^^''''^^^'^^'^W) ^ ikh-<j{t))^r'^f[kh)h. 



h^b 



Remark 2.2. /n Definition ] 2. 8\ we are using summations with limits that are 
reals. For example, the summation that appears in the definition of operator 
a^h^ has the following meaning: 

Y^ G{k) = G{a/h) + G{a/h + 1) + G{a/h + 2) + • • • + G{t/h - v), 

where t G {a + vh, a ^h + vh. a + 2/i + vh. . . . ,a + kh +i'h} with fc G N. 

b 
Lemma 2.9. Let v > Q be an arbitrary positive real number. For any t (zT we 
have: (i) lim.^o a^l" f{t + vh) = f{t); (ii) lim.^o h^b' H^ - ^h) - f{t). 

Proof. Since 

<,A-V(i + ^'^) = f^ [ \t + vh- a{s))\r'^f{s)As 

t_ 

= f^ E (* + ^'^ - <^{kh)ir^^f{kh)h 

ft. — , 



= ^^^^^ + r>TT) E (* + ^'^ - <j{kh)t-^'^f{kh)h , 

h 

it follows that lim,y^o a^.J^'^ f{t + vh) — f{t). The proof of (ii) is similar. D 



For any t e T and for any ly > we define a^lf{t) := hAlJit) := f{t) and 
write 



aA-''f{t + iyh) = h-'f{t) + 



^A- VW = h-'fit - yh) + 



n^+^)Ja 



{t + vh-a{s))^r^^f{s)As, 



<y{b) 



[v-l] 



(8) 



{s + vh~a{t)yr 'fi^)^^ 



Theorem 2.10. Let f G J-j and u > 0. For all t G T" we have 



aA^''f^it + l.h)^UA^''f{t + ,yh)Y 



-{t + iyh-a)\r'^f{a). (9) 



To prove Theorem 12 . f 01 we make use of a technical lemma: 
Lemma 2.11. Let t G T". The following equality holds for all s G T*^; 

A,,,[{t+,.h-s)\r'^f{s))) 

= {t + vh- a{s))i^-'^f^{s) -{v- l){t + vh- a{s)t-^^f{s) 
Proof. Direct calculations give the intended result: 

AsM({t + vh-si;^-^^f{s)) 

= Asj, ((t + iyh- s)i'-'^) f{s) + {t + iyh- ais))i^-'^ f^{s) 



(10) 



fis) 



t-\-uh—a{s) 



r(i±^ + i) 



r 



(^±^^^^^ + i-(.-i)) r(i±f-^ + i -(.-!)) 



+ {t + ^h-o{s)ir'^f''{s) 
h^~^ 






-p^ t-\-iyh — s \ p/ t-\-i'h — s j ^'\ 



n^ + i) Ti^ + 2) 



+ {t + uh-a{s)fr'^f\s) 



n 



t-^vh — s — h 



^^ t-s+i^h-h +i_(^_2)) 



iL-(-(^ - 1)) + {t + :.h- a{s)t-'^ f^{s) 



= -{u -l){t + uh- <j{s)fr'^f{s) + {t + yh - a{s)t-''^ f^{s) , 
where the first equality follows directly from ^. D 

Remark 2.3. Given an arbitrary i G T" it is easy to prove, in a similar way 
as in the proof of Lemma \2.11[ the following equality analogous to (jlOp ; for all 

SGT 



A, 



,,((. + ^/i-a(t))i'^-^V(.s)) 



= {v- l)(s + uh- a{t)fr^^r{s) + {s + vh- ait))i^-'^f^is) . (11) 



Proof, (of Theorem I2.10p From Lemma 12.111 we obtain that 



aA^'^f^it + ^h) = /.V^(i) + Y^;;^ j^ {t + vh~ a{s))t'^ f^{s)As 



h^nt) + 



Tiv + l) 



{u-l) 



{t+vh-s)r"f{s) 



n^+^)Ja 



{u~l){t + uh-a{s)t-^^f{s)As 



(12) 



v{t + vh- a)|;"^^ 



r(^ + i) 



f{a) + h''f''{t) + vh''-^f{t) 



+ 



{i^-l){t + i^h-<j{s))lr^\f{s)As. 



We now show that {aAl^fit + vh))^ equals P^ : 



{,A-^''f[t + yh)Y 



h 



h^!{a{t)) 



r(^ + i)A 



<T(t) 



Ht) + vh-a{s)f^-''^f{s)As 



-h^fit) - 



h'f^it) + 



r(i^ + i)A 



{t + i.h-<7{s))lr'^f{s)As 



hr{i^ + i) 



{t + vh-<j{s))\:-'^f{s)As 



{a{t) + iyh-a{.s))\r'^f{s)As 



h-f^{t) + — ^ f A,,, ((i + ,yh- a(,s))i^-i') /(.)A. + h'^-'i.fit) 



h'f^it) + 



n'^ + ^)Ja 



{u-2) 



{v - l)(i ^vh- a{s)yCf{3)As + vh'^-^fit) 



U 



Fohows the counterpart of Theorem 12.101 for the right fractional h-svnn: 
Theorem 2.12. Let f E Tj and v > 0. For all f e T'' we have 

^A;i;^f^{t-^h) = -^J^^(b + ^h-ait))ir'^fib) + i,A^-^f{t-uh))^. (13) 

Proof. From pil) we obtain from integration by parts (item 2 of Lemma 12. ip 
that 



hA;li)f^{t-,.h) 



r(^ + i) 



f{h) + h'^f'^{t)-vh'^-'f{a{t)) 



r(^ + i)7<.(t) 



{u-l){s + vh^a{t))^;;-^^f%s)As. 



(14) 



We show that {hA'^fit - vh))^ equals (O: 



{uA^'^f{t~uh)Y 



h^rit) + 



hr{i^ + i) 



aib) 



a(b) 



'■(t) 



(i.-l) 



{s + ,yh-a'{t)))t'^fis)As 



a^t) 



h-'f^it) + 



h'^f^it) - 



h'f^it) - 



{s + i'h~ait)Yr'>fis)A. 



vh''-^f{a{t)) 



<T(b) 



<T(b) 



A. 



,,. ((s + vh^ 0(1))^;^-^^) f{s)As - vh''-'f{a(t)) 



r(^ + i)y..(t) 



2/,nn(i'-2) 



(^. - l){s + i./i - a'{t)f;^-'>f{s)As - vh^-^f{cj{t)) 



r(^ + l)^(t) 



(;. - l)(,s + vh- ait))^r'^f{s)As - vh^-^f{a{t)). 



U 



Definition 2.13. Let < a < 1 and set 7 := 1 — a. The left fractional 
difference aAf^f{t) and the right fractional difference hA'^f{t) of order a of a 
function f G jFy are defined as 

aA'f:f{t):^UA--'f{t + jh))^ and ^A^ f{t) := -{hA^yit - ^h))^ 

for allte T«. 



3. Main Results 

Our aim is to introduce the /i-fractional problem of the calculus of variations 
and to prove corresponding necessary optimality conditions. In order to obtain 
an Euler-Lagrange type equation (cf. Theorem 13. 5p we first prove a fractional 
formula of ft,-summation by parts. 

3.1. Fractional h-summation by parts 

A big challenge was to discover a fractional ^.-summation by parts formula 

within the time scale setting. Indeed, there is no clue of what such a formula 

should be. We found it eventually, making use of the following lemma. 

2 
Lemma 3.1. Let f and k be two functions defined on T'^ and T" , respectively, 

and g a function defined on V x T'' . The following equality holds: 



fit) 



git,s)k{.s)As 



•P{b) 

At= I k{t) 



a{t) 



l{s,t)f{s)As 



At. 



Proof. Consider the matrices R = [f{a + h), f{a + 2h), • • • , f{b — h)], 



Ci 



g{a + /i, a)k{a) 
g{a + 2/1, a)k{a) + g{a + 2h,a + h)k{a + h) 

g{b - h, a)k{a) + g{b - h,a + h)k{a + K) A ^ g{b- h,b- 2h)k{b - 2h) 
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Co 



g{a + /i, a) 
g{a + 2h, a) 



, C3 



g{b-h,a) 
Direct calculations show that 

b 



fit) 



g{t,s)k{s)As 





17(0 + 2h, a + h) 

g{b ~ h,a + h) 

b/h-l i-l 



, C4 






g{b-h,b-2h) 



At^h^ Y^ f{ih) Y^ g{ih,jh)k{jh)^h^R-Ci 

i—a/h j—a/h 

h^R ■ [k{a)C2 + k{a + h)C:i + ■ ■ ■ + k{b - 2h)Ci] 

b/h-l b/h-1 

Ha) J2 9iJh,a)fiJh) + kia + h) J^ 9{jh,a + h)f{jh) 

j=a/h+l j=a/h+2 

b/h-1 

+ --- + k{b-2h) Y. 9ijKb-2h)fUh) 

j=b/h~l 
b/h-2 b/h-1 (f,) r fc 

Y H^h)h Y 9ijh,th)f{jh)h^ k{t) / g{s,t)f{s)A^ 

/J, ,.■ /„■^^/I, -J a J<^(t) 



i^a/h 



j=(7(ih)/h 



At. 



n 



Theorem 3.2 (fractional /i-summation by parts). Let f and g be real valued 
functions defined on T'^ and T, respectively. Fix < a < 1 and put 7 := 1 — a. 
Then, 



b np{b) 

f{t)aAlg{t)At = hy{p{b))g{b) - h-'f{a)g{a) + / HA^^^,^f{t)g%t)At 



7 



+ f(7TT)^^''^ (/j* + ^^ " a)(-'-^VWAt - f^^ {t + jh- a{a))i^-'^f{t)At\ 



(15) 



Proof. By © we can write 



/(0aA^5(t)Ai = / fit)UAr9it + lh))^^t 



fit) 



7 



At 



,A-^g-(t + jh) + Yf/^(i + ^^- «)i''"'^5(«) 

I' /(i)aAr5^(i + lh)^t + I' Y^:^^it + jh- a)'-r'^fit)g{a)At. 

(16) 
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Using ^ we get 



b 
fit) 



h'g'^it) + Y^:^^ j^ (t + ^h- <j{s))i^-'^g^{s)As 



At 



h^J^ f{t)g^{t)At + ^,^-^J^ g^{t) J^Js + ^h-a{t))i^-'\f{s)AsAt 



= hy{p{b)Mb) - g{p{b))] + / g^{t)^A;-'f{t - jh)At, 



P(b) 

p(fc)J 

where the third equahty follows by Lemma [01 We proceed to develop the right 
hand side of the last equality as follows: 

hy{p{b)Mb) - g{p{b))] + / .g^(t),A-])/(t - jh)At 
= hy{p{b))[g(b) g{p{b))] + [g{t)nA;^l,^f{t - 7/^)]*^^ ^ 

P(b) 

g^t)UA;^^^f{t~jh))^At 



hy{p{b))g{b) - hy{a)g{a) 

M I [s + lh- a{a))^^-'^f{s)As + f („A^(,)/(t)) g'^{t)At, 



7 



r(7+l)"' ' J.(a) 

where the first equality follows from Lemma l2.ll Putting this into (1161) we get 
(Ha. D 



3.2. Necessary optimality conditions 

We begin to fix two arbitrary real numbers a and /3 such that a, /3 e (0, 1]. 
Further, we put "f :— 1 — a and v := 1 — (3. 

Let a function L{t,u,v,w) iT^xMxMxM— j^Mbe given. We consider 
the problem of minimizing (or maximizing) a functional C : Ff — J> M subject to 
given boundary conditions: 

^y{-)) = / Ht, y%t),aA'-^y{t),hAly{t))At -^ min, y{a) = A, y{b) - B . 

J a 

(17) 
Our main aim is to derive necessary optimality conditions for problem ([T7)) . 

Definition 3.3. For f ^ Fj we define the norm 

||/|Hmax|r(i)|+maxUA^/(i)|+maxUAf/(t)|. 

A function y G J-f with y{a) — A and y{b) = B is called a local minimum, for 
problem (|17p provided there exists (5 > such that C{y) < C{y) for all y G J-j 
with y{a) — A and y{b) — B and \\y — y\\ < 5. 

Definition 3.4. A function rj G J-j is called an admissible variation provided 
r] ^ and r]{a) = 77(6) = 0. 

From now on we assume that the second-order partial derivatives £„„, L^v, 
Luw, Lyj^, Lyy, and i^tu exist and are continuous. 
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3.2.1. First order optimality condition 

Next theorem gives a first order necessary condition for problem (|17l) . i.e., 
an Euler-Lagrangc type equation for the fractional /i-difFerence setting. 

Theorem 3.5 (The /i- fractional Euler-Lagrange equation for problem pT)) ). // 
y G J-'j is a local minimum for problem (|17p . then the equality 

Lu[y]{t) + hA'^pibMm) + aKL^m) = (18) 

holds for allt G T'^ with operator [■] defined by [y]{s) = (s, y'^(s), QA"y(s), sAjjy(s)). 

Proof. Suppose that y{-) is a local minimum of C[-]. Let ri{-) be an arbitrarily 
fixed admissible variation and define a function $ : ( — j-ryiT , yT)\\ ) ~^ I^ by 

^e) = C[y{-) + evi-)]. (19) 

This function has a minimum at £ = 0, so we must have $'(0) = 0, i.e., 



Lu[m)v''it) + L,[y]{t)aA"^v{t) + L^[m)h^Ht) 



At = 0, 



which we may write equivalently as 

rp(b) 



hLu[mW{t)\t^p[b)+ L4y]it)r,''it)At+ LM{t)aK^(t)M 

J a J a 

+ / L^mt)hAU{t)At^O. (20) 

J a 

Using Theorem 13.21 and the fact that 77(a) = ry(6) = 0, we get 

b rp{b) 

L4m)aKvit)At ^ J (;,A^(,)(L4y])(i))rr(t)Ai (21) 

for the third term in (j20p . Using (J13p it follows that 



L4y]{t)hA^^r,{t)At 



L^[y]{t){hA-''ri{t-i,h))^At 



Ln,[ym 



At 



r L4ym>^^;n)V^it - ^h)At + ,0^^ j\h + uh~ a{t))t~'''L^[y]{tW. 

(22) 
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We now use Lemma 13.11 to get 

l-b 

L^[m)hA-i;^rj'^{t-i,h)At 



J a 



h'^v'^it) + 



r(^ + i)^(t) 



{s + uh-a{t))^r^\^{s)As 



At 



/ /i^L^[y](t)77^(<)At 

J a 



n'^+i)Ja 



Ja(t) 



At 



h'L^ymv'^WAt 



+ 



r(r^ + i)7a 



b r 



(..-1). 



r^^it) / {t + ,.h-ai.s)yr >L^mis)A 



At 



Tj'^{t)aA^''{L4y]){t + j,h)At. 



(23) 



We apply again the time scale integration by parts formula (Lemma 12. ip , this 
time to (P5|) . to obtain, 



7y^(i)aA^''(L^,[y])(f + z//i)Ai 

•p(&) 

V^it)aA-''iL^[y]){t + :yh)At 

+ (77(6) - fjipibjJUA--- [L^y]) {t + ^/i)|,=p(b) 



= [ry(t),A^''(L^[y])(t + ^/i)];:f ^ - / rj^t)UAj;'' {L^[y]) {t + ,.h))^At 

J a 

+ v{b)aA-'' {L^[y]) {t + i^h)\t^p^h) - v{p{b))aA--' {L^[y]) {t + Wi)|t^p(fc) 
= vib)aA^'' {L^m (t + ;^/i)|t=p(fc) - via)aA^-' {L^[y]) [t + vh)\t=a 

^'^{t)aAl{L^[y]){t)At. 

(24) 



[t)At, 



Since 77(a) = 77(6) = we obtain, from ([23)1 and ((24)) . that 
/ L^[y](f),A;(^)77^(t)Af = - / 7?-(t),Af(L, 

-'a J a 

and after inserting in ([2^ . that 

L^[y](i)„Af77(t)Ai= / 7;^(t),Af (L^[y]) (i)At. 
By (ini) and (gS]) we may write (HO]) as 

U[y\{t) + ,A^(,) (L„[y]) [t) + aAf (L„[y]) (i) 77-(t)Ai = . 



(25) 
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Since the values of rj"' (i) are arbitrary for t E T'^ , the Euler-Lagrange equation 
dH]) holds along y. D 

The next result is a direct corollary of Theorem 



Corollary 3.1 (The /i-Euler-Lagrange equation - cf., e.g., 18ll24|)- LetT be the 



time scale /iZ, h > 0, with the forward jump operator a and the delta derivative 
A. Assume a,b Cz T, a < b. If y is a solution to the problem 



b 

A 



C{y{-))^ / Lit,y-{t),y^it))At^m[n, 2/(a) = A, y{b) = B , 



then the equality Lu{t,y'^{t),y'^{t)) - {L^{t,y'^{t),y'^{t))) = holds for all 

ter\ 

Proof Choose a ~ 1 and a L that does not depend on w in Theorem 13.51 D 
Remark 3.1. // we take h = I in Corollary \3.1\ we have that 
Luit, r{t), Ay{t)) - Ai.(i, r{t),Ay{t)) = 

holds for all t E T'^ . This equation is usually called the discrete Euler-Lagrange 
equation, and can be found, e.g., in \2!a . Chap. 8]. 

3.2.2. Natural boundary conditions 

If the initial condition y{a) = ^ is not present in problem (J17p (i.e., y{a) 
is free), besides the /i-fractional Euler-Lagrange equation ([TSt the following 
supplementary condition must be fulfilled: 

{t + jh-a{a))[^-^^L4y]{t)^t] +^»[y](a) = 0. (26) 

<T(a) / 

Similarly, if y{b) = B is not present in (|17l) {y{b) is free), the extra condition 
hL^[y]{p{b)) + h-^LMipih)) - h^L^[y]{p{b)) 

b 



ib + iyh-a{t))i'-'^L4y]{t)At 



{p{b) + iyh-a{t))[^~'^L4y](t)^tj =0 (27) 



is added to Theorem 13.51 We leave the proof of the natural boundary conditions 
(PB)l and (|77|) to the reader. We just note here that the first term in (l?7|) arises 
from the first term of the left hand side of (PO]) . 



15 



3.2.3. Second order optimality condition 

We now obtain a second order necessary condition for problem ([T7|. i.e., 
we prove a Legendre optimality type condition for the fractional ft.-difFerence 
setting. 

Theorem 3.6 (The ft,- fractional Legendre necessary condition). If y ^ J-t is a 
local minimum for problem (|17[) . then the inequality 

h^L^uMt) + 2/i^+iL™[y](i) + 2ft'^+i(i/ - l)Lu^mt) + h^^il - lfL,M{^{t)) 
+ 2ft'^+^(7 - l)L,^[y]{a{t)) + 2h''+^{u - l)L,^[y]{t) + h^'i^ - lfL^w[m) 

As 



h''L^,[y]{t)+ f h^L,M{s){^r^^{s + ^h~a{a{t)))'^;:-A As > 

Ja{a{t)) VU7 + lj / 

(28) 



holds for allte T" ,, where [y]{t) = {t,rit), aAfy{t),tA^y{t)). 

Proof. By the hypothesis of the theorem, and letting $ be as in (J19p , we have 
as necessary optimality condition that $"(0) > for an arbitrary admissible 
variation ?/(•). Inequality $"(0) > is equivalent to 

b , 

LuuimtWit)? + 2LuM{tK{t)a^lv{t) + 2Lu^[y]{tK{t)hAlr^{t) 



+L,M{t){aAlT^{t)f + 2L,^[y]{t)aA'^vit)hA^bVit) + L„^{t){hAflii{t)f 



At > 0. 
(29) 



Let r e T*^ be arbitrary, and choose 77 : T ^- M given by 77(t) = < „ ,, . 

■^ ' ^ -^ '^ ^ [ otherwise. 

It follows that 77(a) = 77(6) ~ 0, i.e., 77 is an admissible variation. Using (jH]) we 

get 

[Luu[mW{t)f + 2UM{tK{t)aKv{t) + i™[y] W(aA^77(i))2] At 
rb ' 



Luuimwit))' 

+ 2L^4y]it)ri''it) (^h''rj^{t) + ^^^^^ j\t + jh- ais))i^-'\^{s)As^ 
+ L,M{t) (^ftVW + j;^:^^*(< + 7/»-a(s))(''"V(s)A.s) At 

+ £^ L.„[y](t)(^/7V(i) + j;(:^/^'(t + 7/^-^W)i""'^^^(.'^)A.s) At. 
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Observe that 



r-fc 

Lvv[v]{t) 



aHr) 



7 



r(7 + i)ya 

7 



(t + 7;i-a(5))J7"'^77^(s)As) At 



^'^ i™[y](t)(/iVW + ^T^:^£(i + 7/i-<TW)i-^-^)r7^(s)As) At. 

Let t e [a^{T),p{b)] n /iZ. Since 



JT(:^ y^ (^ + 7/^ - ^(.))i-^-'^r;^(.)A. 



7 



r(7 + i) 



/i:;^7 r 

r(7 + i) 
7/1^ 
r(7 + i) 



<T(r) 



(7-l)^A/ 



+ / (t + 7/i-T(5))|7-')77^(s)As 
{t + jh- cT{r))'^r'^ -{t + ^h a{a{T))ir'\ 

(Sr+7-i)r(^+7-i)-(Sr)r(Sr + 7-i) 
(^)r(^) 



-^^f|^(^ + 7/.-.(.(r))r^ 



(30) 



we conclude that 






At. 



Note that we can write t^lvi^) ~ ^h^ (t)'?^!^ ^ ^^) because 77(6) = 0. It is 
not difficult to see that the following equality holds: 



2i„^[y](t)r7"W/iA^r;(t)At 



2Lun.[y]{tK{t)h^7b)'n''i^ - ^h)^t 



p{b)' 



= 2h^+''L„ 



-)(^-l). 



Moreover, 

rb 



2L^^[y]{t)aA'^T]{t)hA^T]{t)At 



2 / L,^[m)\ {h^V^{t) + 



r(7 + l) Ja 



(t + jh-<j(s))^^-'\^{s)As 



/i''r;^(t) + 



{s + iyh^a{t))[^-^\^is)As 



2h''+''+\u - l)L,^[y]{T) + 2/i^+''+i(7 - l)i,^[y](a(T)). 



At 
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Finally, we have that 

rb 



Lww[y]it) 



h-^rj^it) 



n^+i)j.it) 



{s + i^h-a{t))\'^-'^r]^{s)As 



At 



Lww[y]{t) 



-\ 2 






At 



+ hL^^mirW - ^h^Y + h'''+'L^^[y]{a{r)) 



7^ {(.r + ,yh- ait))t'^ - {a{T) + i^h - a{t))\r'^ } 






Similarly as we did in ((30)) . we can prove that 



^W+l) {(^ + ^^ - '^W)^'^ - ('^(^) +^h- cj{t)t-'^ } 



h^'^^ir +.h-mfr\ 



Thus, we have that inequality (j29l) is equivalent to 



hi h^LuuMt) + ^h-'+^LuMit) + h-'Lv^mt) + L™('T(t))(7/i^ - /iT)2 






As 



+ 2h''+^L^M{t){u - 1) + 2;j^+'^(t. - l)L,^[y](t) 

+ £ /i^L^^[y](s) (^|^^(i + ^h- 'yi.s))^-^^ aX > 0. (31) 

Because h > 0, (|5T|) is equivalent to (PSJ) . The theorem is proved. □ 

The next result is a simple corollary of Theorem 13.61 

Corollary 3.2 (The /i-Legendre necessary condition - cf. Result 1.3 of 18^). 
Let T be the time scale KL, h > 0, with the forward jump operator a and the 
delta derivative A. Assume a, 6 e T, a < b. If y is a solution to the problem 

nb 

C{y{-))^ / L{t, y^t),y'^{t)) At ^min, y{a) ^ A, y{b) ^ B , 

J a 

then the inequality 

h^L^uMt) + 2hLuv[y]{t) + LvvivKt) + Lvv[y]icT{t)) >0 (32) 

holds for all t e T'^', where [y](t) = {t,y''{t),y^{t)). 
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Proof. Choose a — 1 and a Lagrangian L that does not depend on w. Then, 
7 = and the result fohows immediately from Theorem 13.61 D 



Remark 3.2. When h goes to zero we have cr(t) = t and inequality (1321) coin- 
cides with Legendre's classical necessary optimality condition L^y[y]{t) > (cf., 
e.g., fM 



4. Examples 

In this section we present some illustrative examples. 
Example 4.1. Let us consider the following problem: 

Civ) - \ 



oAlyit)) Af^min, y(0) = , 2/(1) = 1 . 



(33) 



We consider i33\) with different values of h. Numerical results show that when 
h tends to zero the h-fractional Euler- Lagrange extremal tends to the fractional 
continuous extremal: when h ^ i33\) tends to the fractional continuous vari- 
ational problem in the Riemann-Liouville sense studied in jj, Example 1], with 
solution given by 

. N 1 /" dx , , 

y^^^ = 9 / 7. — 7 — ^ ■ (^4) 

^ Jo [(l--x){t-x)]^ 
This is illustrated in Figure [H Ln this example for each value of h there is 




Figure 1: Extremal y{t) for problem of Example 14. II with different values of h: h = 0.50 (•); 
h = 0.125 (+); h = 0.0625 (*); h = 1/30 (x). The continuous line represent function 13411 . 



a unique h-fractional Euler- Lagrange extremal, solution of (|18p . which always 
verifies the h-fractional Legendre necessary condition 

Example 4.2. Let us consider the following problem: 

r-l 



^y) 



i(oA^y(t))2-y-(t) 



At 



2/(0) = 0, 2/(1) =0. (35) 



We begin by considering problem i35\) with a fixed value for a and different 
values of h. The extremals y are obtained using our Euler- Lagrange equation 
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U8\). As in Example \4- l\ the numerical results show that when h tends to zero the 
extremal of the problem tends to the extremal of the corresponding continuous 
fractional problem of the calculus of variations in the Riemann-Liouville sense. 
More precisely, when h approximates zero problem \3^) tends to the fractional 
continuous problem studied in l£, Example 2]. For a ~ \ and h ^ the extremal 
of !i35]) is given by y{t) = \t{l — t), which coincides with the extremal of the 
classical problem of the calculus of variations 



c{y) 



^\y'{tf-yit)]dt^min, y{0) = , y{l) = . 



This is illustrated in Figure\^ for h — -^ , i ^ 1,2,3,4. In this example, for 




Figure 2: Extremal y{t) for problem II35I I 
with a = 1 and different values of h: h = 
0.5 (•); h = 0.25 (x); h = 0.125 (+); h = 
0.0625 {*). 



Figure 3: Extremal y{t) for II35I I with h = 
0.05 and different values of a: a = 0.70 
(•); a = 0.75 (x); a = 0.95 (+); a = 0.99 
(*). The continuous line is j/(t) = 2^(1—*). 



each value of a and h, we only have one extremal (we only have one solution 
to il8\) for each a and h). Our Legendre condition (I28p is always verified along 
the extremals. Figure\^ shows the extremals of problem psp for a fixed value of 
h (h = 1/20J and different values of a. The numerical results show that when 
a tends to one the extremal tends to the solution of the classical (integer order) 
discrete-time problem. 

Our last example shows that the /i-fractional Legendre necessary optimality 
condition can be a very useful tool. In Example 14.31 we consider a problem 
for which the /i-fractional Euler-Lagrange equation gives several candidates but 
just a few of them verify the Legendre condition ((28)) . 



Example 4.3. Let us consider the following problem: 



^y) 



iaA^yit)f + 9UAtyit)fAt 



y(a) = 0, y(fc)-l 



(36) 

For a = 0.8, /3 = 0.5, h = 0.25, a = 0, 5 = 1, and 9 = 1, problem ^ has 
eight different Euler-Lagrange extremals. As we can see on Ta&/e[7] only two of 
the candidates verify the Legendre condition. To determine the best candidate 
we compare the values of the functional C along the two good candidates. The 
extremal we are looking for is given by the candidate number five on Table QJ 
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# 


vii) 


va) 


Vii) 


C{y) 


Legendre condition II28II 


1 


-0.5511786 


0.0515282 


0.5133134 


9.3035911 


Not verified 


2 


0.2669091 


0.4878808 


0.7151924 


2.0084203 


Verified 


3 


-2.6745703 


0.5599360 


-2.6730125 


698.4443232 


Not verified 


4 


0.5789976 


1.0701515 


0.1840377 


12.5174960 


Not verified 


5 


1.0306820 


1.8920322 


2.7429222 


-32.7189756 


Verified 


6 


0.5087946 


-0.1861431 


0.4489196 


10.6730959 


Not verified 


7 


4.0583690 


-1.0299054 


-5.0030989 


2451.7637948 


Not verified 


8 


-1.7436106 


-3.1898449 


-0.8850511 


238.6120299 


Not verified 



Table 1: There exist 8 Euler-Lagrange extremals for problem I I36I I with a = 0.8, /3 = 0.5, 
h = 0.25, a = 0, 6 = 1, and 6 = 1, but only 2 of them satisfy the fractional Legendre condition 



# 


5(0.1) 


5(0.2) 


5(0.3) 


5(0.4) 


^(5) 


IMli 


1 


-0.305570704 


-0.428093486 


0.223708338 


0.480549114 


12.25396166 


No 


2 


-0.427934654 


-0.599520948 


0.313290997 


-0.661831134 


156.2317667 


No 


3 


0.284152257 


-0.227595659 


0.318847274 


0.531827387 


8.669645848 


No 


4 


-0.277642565 


0.222381632 


0.386666793 


0.555841555 


6.993518478 


No 


5 


0.387074742 


-0.310032839 


0.434336603 


-0.482903047 


110.7912605 


No 


6 


0.259846344 


0.364035314 


0.463222456 


0.597907505 


5.104389191 


Yes 


7 


-0.375094681 


0.300437245 


0.522386246 


-0.419053781 


93.95316858 


No 


8 


0.343327771 


0.480989769 


0.61204299 


-0.280908953 


69.23497954 


No 


9 


0.297792192 


0.417196073 


-0.218013689 


0.460556635 


14.12227593 


No 


10 


0.41283304 


0.578364133 


-0.302235104 


-0.649232892 


157.8272685 


No 


11 


-0.321401682 


0.257431098 


-0.360644857 


0.400971272 


19.87468886 


No 


12 


0.330157414 


-0.264444122 


-0.459803086 


0.368850105 


24.84475504 


No 


13 


-0.459640837 


0.368155651 


-0.515763025 


-0.860276767 


224.9964788 


No 


14 


-0.359429958 


-0.50354835 


-0.640748011 


0.294083676 


34.43515839 


No 


15 


0.477760586 


-0.382668914 


-0.66536683 


-0.956478654 


263.3075289 


No 


16 


-0.541587541 


-0.758744525 


-0.965476394 


-1.246195157 


392.9592508 


No 



Table 2: There exist 16 Euler-Lagrange extremals for problem II36I I with a = 0.3, h = 0.1, 
a = 0, b = 0.5, and 9 = 0, but only 1 (candidate #6) satisfy the fractional Legendre condition 
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For problem h36\) with a = 0.3, h = 0.1, a ~ Q, b = 0.5, and 6 — 0, we 
obtain the results of Table O' there exist sixteen Euler- Lagrange extremals but 
only one satisfy the fractional Legendre condition. The extremal we are looking 
for is given by the candidate number six on Table [H 

The numerical results show that the solutions to our discrete-time fractional 
variational problems converge to the classical discrete-time solutions when the 
fractional order of the discrete-derivatives tend to integer values, and to the 
fractional Riemann-Liouville continuous-time solutions when h tends to zero. 

5. Conclusion 

The discrete fractional calculus is a recent subject under strong current de- 
velopment due to its importance as a modeling tool of real phenomena. In this 
work we introduce a new fractional difference variational calculus in the time- 
scale {hZ)a, h > and a a real number, for Lagrangians depending on left and 
right discrete-time fractional derivatives. Our objective was to introduce the 
concept of left and right fractional sum/difference (cf. Definition 12. 8p and to 
develop the theory of fractional difference calculus. An Euler-Lagrange type 
equation ([T5)) . fractional natural boundary conditions (J26D and (P7|) . and a sec- 
ond order Legendre type necessary optimality condition (|28p. were obtained. 
The results are based on a new discrete fractional summation by parts formula 
(|15p for {hZ)a. Obtained first and second order necessary optimality conditions 
were implemented computationally in the computer algebra systems Maple and 
Maxima. Our numerical results show that: 

1. the solutions of our fractional problems converge to the classical discrete- 
time solutions in {h'L)a when the fractional order of the discrete-derivatives 
tend to integer values; 

2. the solutions of the considered fractional problems converge to the frac- 
tional Riemann-Liouville continuous solutions when /i ^^ 0; 

3. there are cases for which the fractional Euler-Lagrange equation give only 
one candidate that does not verify the obtained Legendre condition (so 
the problem at hands does not have a minimum); 

4. there are cases for which the Euler-Lagrange equation give only one can- 
didate that verify the Legendre condition (so the extremal is a candidate 
for minimizer, not for maximizer) ; 

5. there are cases for which the Euler-Lagrange equation give us several 
candidates and just a few of them verify the Legendre condition. 

We can say that the obtained Legendre condition can be a very practical tool to 
conclude when a candidate identified via the Euler-Lagrange equation is really 
a solution of the fractional variational problem. It is worth to mention that a 
fractional Legendre condition for the continuous fractional variational calculus 
is still an open question. 

Undoubtedly, much remains to be done in the development of the theory 
of discrete fractional calculus of variations in {h'Z)a here initiated. Moreover, 
we trust that the present work will initiate research not only in the area of 
the discrete-time fractional calculus of variations but also in solving fractional 
difference equations containing left and right fractional differences. One of the 
subjects that deserves special attention is the question of existence of solutions 
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to the discrete fractional Euler-Lagrange equations. Note that the obtained 
fractional equation (IT5| involves both the left and the right discrete fractional 
derivatives. Other interesting directions of research consist to study optimalitj' 
conditions for more general variable endpoint variational problems 28|, [SJ, 



isopcrimctric problems fS, 65 ; higher-order problems of the calculus of variations 
[12, 24, 38] ; to obtain fractional sufficient optimality conditions of Jacobi type 
and a version of Noether's theorem [17|, |2l|, |25|, |27[ for discrete-time fractional 
variational problems; direct methods of optimization for absolute extrema fl_9|, 
36l . l49j ; to generalize our fractional first and second order optimality conditions 
for a fractional Lagrangian possessing delay terms [IJ, |37| ; and to generalize the 
results from {h'Z)a to an arbitrary time scale T. 
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